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Abstract 

The various roles of boundary terms in the gravitational Lagrangian and 
Hamiltonian are explored. A symplectic Hamiltonian-boundary-term ap- 
proach is ideally suited for a large class of quasilocal energy-momentum ex- 
pressions for general relativity. This approach provides a physical interpreta- 
tion for many of the well-known gravitational energy-momentum expressions 
including all of the pseudotensors, associating each with unique boundary 
conditions. From this perspective we find that the pseudotensors of Einstein 
and M0ller (which is closely related to Komar's superpotential) are especially 
natural, but the latter has certain shortcomings. Among the infinite possi- 
bilities, we found that there are really only two Hamiltonian-boundary-term 
quasilocal expressions which correspond to covariant boundary conditions; 
they are respectively of the Dirichlet or Neumann type. Our Dirichlet expres- 
sion coincides with the expression recently obtained by Katz and coworkers 
using Noether arguments and a fixed background. A modification of their 
argument yields our Neumann expression. 
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I. INTRODUCTION 



Via their energy-momentum density, material sources generate gravitational fields. 
Sources interact with the gravitational field locally, hence they should be able exchange 
energy-momentum with the gravitational field — locally. From this physical conception 
we are led to expect the existence of a local density for gravitational energy-momentum. 
However the identification of a good localization for gravitational energy-momentum has 
turned out to be an outstanding fundamental problem. Standard techniques led only to 
various reference frame dependent complexes referred to as pseudotensors. This result can 
be understood in terms of the equivalence principle which implies that one can not detect 
any feature of the gravitational field at a point. Consequently, the whole idea has been 
criticized (see, e.g., [^], p 467) and the pseudotensor approach in particular has largely 
been abandoned. 

A new idea: quasilocal (i.e., associated with a closed 2-surface) has become widely 
accepted. Many quasilocal proposals have been considered (for the older works see [§; 
more recent works are cited in |ll4| , p!0| ). Many well-known quasilocal expressions obtained 
by different approaches have been discussed in the literature; although they generally give 
different values , most seem to work well enough at least for certain physical situations. A 
number of criteria for selecting a good quasilocal expression (see, e.g., |TH]), including good 



limits at spatial infinity, at future null infinity, to weak fields, and to fiat spacetime have 
been advocated. Such requirements, however, have proved to be insufficient; in fact it has 
been noted that there still exist an infinite number of expressions satisfying the proposed 
criteria f^. We infer that additional principles and criteria are very much needed to reduce 
and to parameterize, if not entirely eliminate, the ambiguity. 

One might hope that there would exist a "best" gravitational energy-momentum ex- 
pression which has either not yet been identified or at least not yet accorded widespread 
acceptance. On the other hand, it is well to keep in mind that there are physical situa- 
tions where there is not one unique energy. One example is thermodynamics, wherein there 
are several energies (viz., internal, enthalpy, Gibbs and Helmholtz) corresponding to differ- 
ent choices of boundary conditions and independent variables; each one gives the relevant 
value of the energy for a particular physical situation. An even more appropriate exam- 
ple is electrostatics. It is well known that the work done in moving a system of charges 
and dielectrics differs depending on whether one holds the potential or the charge density 
fixed. We expect gravity to behave in a similar fashion: consequently various definitions 
of gravitational energy may each he associated with their own unique boundary condition. 
Fortunately, there exists a systematic technique, symplectic analysis |^ , which, along with 
its associated control-response relations, can be used to identify the relationship between an 
energy-momentum expression and the associated boundary conditions. 

Our approach to quasilocal energy-momentum is by way of the Hamiltonian formulation 
— essentially, we take energy to be given by the value of the Hamiltonian. The rationale goes 
back to Noether's work connecting symmetries and conserved currents: in particular energy- 
momentum is associated with translations in spacetime. The Hamiltonian is the Noether 
canonical generator of timelike displacements. The Hamiltonian for gravitating systems 
for a finite region of spacetime includes, in addition to a volume density term, a surface 
term which plays a key role. Its value will determine the quasilocal energy-momentum and, 
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through its variation, it governs the associated boundary conditions. In an earher work 



IJ] we presented our ideas, apphed to rather general geometric gravity theories including 
Einstein's general relativity (GR), in terms of differential forms. However that technique is 
not so well known by many people interested in gravitational energy-momentum; moreover 
most people are mainly interested specifically in GR. Hence we present here a formulation 
of the application of our ideas to GR in the more traditional holonomic (coordinate basis) 
style. This will serve not only to make our ideas more widely accessible but will also facilitate 
comparison with the results obtained by other investigators, e.g., []7|, p!8| , p6| ,|8|, |43|j39[| . 

Elsewhere we have used our Hamiltonian-boundary-term approach to quasilocal quan- 
tities to show that all of the pseudotensors give well defined quasilocal energy-momentum 



values, each of which is associated with a particular choice of boundary conditions ||T^ . From 
the Hamiltonian boundary approach we find that two, the Einstein and M0ller expressions, 
each arise naturally, although the latter has a few extra shortcomings. A more important 
failing however, in our opinion, is that none of the pseudotensors is associated with a truly 
covariant boundary condition. 

For the required new principle and criteria to restrict the GR quasilocal energy- 
momentum expression we have advocated the Hamiltonian boundary variation principle and 
the criteria of covariance. For GR we found that there are only two covariant quasilocal ex- 
pressions (both of which depend on the choice of a reference configuration on the boundary) 
which correspond, respectively, to Dirichlet and Neumann boundary conditions. At first we 
were surprised to learn that our Dirichlet expression coincides with an expression developed 



by Katz and coworkers ||3^ , |28| , P9| using a different approach based on a Noether type argu- 
ment at the Lagrangian level along with a fixed global background geometry. With hindsight 
we now see that this agreement is related to the close connection between the Hamiltonian 
and Noether translation current. Here we show how to modify their argument to also obtain 
our Neumann quasilocal energy-momentum expression. Along the way we clarify the roles 
of the boundary (or total derivative) terms in the Lagrangian, the Hamiltonian and their 
respective variations. 

II. THE SYMPLECTIC IDEA IN GENERAL RELATIVITY 

In this section we outline the symplectic idea for Lagrangian and Hamiltonian formula- 
tions (here and elsewhere we are much influenced by Kijowski and coworkers [^,^^-^) for 
general relativity (GR), Einstein's theory of gravity (a detailed discussion for general geo- 
metric gravity theories, in terms of differential forms, appears in [Q). The simple and direct 



way to reveal the symplectic structure of a physical configuration is through the variation 
of the associated Lagrangian or Hamiltonian. 

A. Lagrangian formulation 

Let us first briefly review some features of the Lagrangian variational principle for clas- 
sical field theories. For our purposes we found it convenient to use the first order formalism. 
For a field (j)^ the first order Lagrangian scalar density has the form 

c = P!^d,<p^-A{<p^,P!^). (1) 
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The field equations (whicli in this formulation contain only first derivatives of the fields) are 
taken to be the Euler-Lagrange expressions implicitly determined by the variation: 



5C = d,{P!t6<p^) + j-5<p^ + j^5Pl (2) 



The action is obtained by integrating the Lagrangian density over a spacetime region; the 
variation of the action is given by the integral of (0). When integrated over a spacetime 
region the total derivative term becomes a boundary term. Technically the variational 
derivatives of the action are well defined only if the boundary term in the variation vanishes. 
That requirement shows what must necessarily be held fixed on the boundary, this quantity 
is referred to as the "control variable". In this case it is the field 0"^, thus this Lagrangian is 
differentiable only on the space of fields for which is given a predetermined dependence on 
the boundary. The variation boundary term, moreover, has a certain symplectic form which 
connects the "control variable" with an associated "response variable" (in this case -Pji). 

We are only concerned with actions which do not depend on the position except via the 
fields. Hence they have an invariance under local infinitesimal translations (i.e, diffeomor- 
phisms) — which can be represented by Lie derivatives. Thus, for an arbitrary vector field 
AT, 



6C ^ 6jC_ 



£mC := d.iN'^C) = d^{P>X£Mr) + T-A^Nr + -t1^£nP^a- (3) 

od)^ or. 



From this identity we conclude, by taking A^'^ to have constant coefficients, that the canonical 
energy-momentum density, 

:= P>id,ct>^ - b^,L, (4) 

is conserved. More precisely its divergence is proportional to a combination of the field equa- 
tions and hence vanishes "on shell". Note that the canonical energy-momentum density is 
not unique in the sense that we can add to it an expression which is automatically diver- 
gence free. Such an expression is necessarily of the form d^Uv^"' where Uy^"^ = —Uy"'^. This 
ambiguity allows one to adjust the zero of energy and has been exploited to find "improved" 
energy-momentum tensors such as the symmetrized one constructed by Belinfante and 
Rosenfeld ETl . On the other hand, since we have also assumed invariance with non-constant 



A^'^, we are also requiring that @ be identically satisfied for the terms proportional to dN 
(this is only possible if the list of dynamic variables includes certain geometric variables). 
In this way we discover that T^j^ itself is linear in the field equations and thus vanishes "on 
shell" . In other words the 'conservation law' is actually a differential identity connecting the 
field equations showing that they are not all independent, hence the evolution of the dynam- 
ical variables is underdetermined — a fact which is directly related to the local 'translational' 
gauge (i.e., diffeomorphism) freedom of the theory. 

Now let us apply this analysis to gravity. There are several choices of variables and 
Lagrangians which can be used. Since we favor a first order approach, a natural geomet- 
ric choice is to regard the metric and connection as independent fields. Even within this 
overall approach there are various options, in particular the metric degrees of freedom can 
alternately be encoded in terms of an orthonormal frame while the torsion free and metric 
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compatibility conditions can be imposed a priori, or enforced via Lagrange multipliers, or 
they can be obtained as dynamic field equations (this is easily done in the vacuum case 
which is all that we consider here). All of these approaches merit consideration. We have 
investigated many of the possible combinations; our preliminary conclusion is that they lead 
to essentially the same result ||^. Here we consider explicitly only one case which is relatively 
simple in the holonomic treatment. 

The field equations of (vacuum) GR can be obtained from a first order variational prin- 



ciple using the Hilbert Lagrangian density in the so-called 'Palatini' form |^| 



CH:=Tr^''RUr). (5) 

Here we are using the contravariant metric density, defined by tt^'^ := {2n)^^y/—gg^^'^ , where 
K := SttG/c^, and the (symmetric, i.e., torsion free) connection coefficients, T^ctf3, as inde- 
pendent variables (while our conventions are generally those of |^^, our treatment can be 
compared with |30|j3T| in which the same variable combinations appear). The variation of 
the Lagrangian density, after the usual integration by parts, has the form 

5Ch = ^^TT + + {7:^^51: 5T\p) . (6) 

The variational derivatives will give the desired field equations: R^v = 0, from the variation 
with respect to n^'^, and Dxh'^^ = (equivalent to Dxg^^ = 0), from the variation with 
respect to F. When integrated over a spacetime region, the total derivative term gives rise 
to a boundary term. This boundary term shows that the control variable is the connection 
and the response variable is linear in tt^'^. The variational derivatives are well defined only if 
this boundary variation term vanishes For a finite region this means we must 'control' or 



'hold fixed' (i.e., give as a prespecified function) F on the boundary. For an asymptotically 
fiat region the connection vanishes asymptotically, nevertheless F = is not a sufficient 
boundary condition. Since we must allow for variations with the generic spatial fall offs 
671 ~ 0(l/r), (5F ~ 0(l/r^), the Lagrangian boundary variation term will yield a finite 
result in the asymptotic limit. Formally the situation is then described by saying that, in 
this case, the variational derivatives of the action are not well defined on the full space 
of asymptotically flat metric and connection fields, but rather only on the subspace where 
we actually fix the specific asymptotic form of F. This 'problem' is closely related to the 
fact that the Hilbert Lagrangian density is asymptotically 0(1/ r^); consequently the action 
diverges for r 00. The remedy is simple: adjust the Lagrangian density by adding a total 
derivative term. 

For GR, an obvious alternative is the "first order" (in derivatives of the metric) La- 
grangian density, initially introduced by Einstein, which can be easy obtained by adding a 
total derivative term to (ISI): 



Ce := Ch + d, (vr^'^F^^^^-j) 

= TT^-^ (f^^F'^,^ - F%,F^„^) . (7) 



The variation of the Einstein Lagrangian, 

5Ce = (fields eq. terms) + (f^„^5;^J dn^") , (8) 
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has the same field equation terms as before but a different boundary term which reflects an 
alternate symplectic structure and shows that the variable to be held flxed — the "control 
variable" — is now the contravariant metric density while the "response variable" is a cer- 
tain combination of the connection. Now, for asymptotically flat fall offs, the Lagrangian 
boundary variation term does vanish asymptotically, so the variational derivatives are well 
deflned on the space of all asymptotically flat flelds (a related fact is that the Einstein action 
is flnite). 

However the big drawback is that we now have a response expression which is linear in F, 
a non-tensorial, reference frame dependent object (along with this the Lagrangian density 
itself is not covariant). The cure for this new 'problem' is to introduce a background (or 
reference) connection f (actually this is really only essential on the boundary) and deflne 
AF := F — F. The latter, being the difference between two connections, is a covariant 
quantity, which can be used in the Lagrangian density boundary (i.e., total derivative) 
term. The 'improved Einstein' action is now 

CiE = CH + d, (vr^'^AF^^^-j) . (9) 

The variation gives the same fleld equation terms but now has a covariant boundary- variation 
symplectic structure: 

5CiE = (flelds eq. terms) + (AF^„^5°J Sti^'') . (10) 

The Lagrangian boundary term does not, as is well known, affect the fleld equations. 
What it does affect is the boundary conditions implicit in the action. From another point of 
view, changing the action by a total derivative term amounts to a canonical transformation; 
in particular, as we saw in the cases considered, it is possible to interchange the role of 
'coordinate' and 'momentum'. Thus the Lagrangian variational boundary term possesses 
important information: the symplectic structure representing the control-response relation 



of the system [32,45|. For instance, the symplectic structure in {m shows that the connection 



is the control variable and the response is a certain combination of the metric. 



B. Hamiltonian formulation 

The energy of a gravitating system can be identifled with the value of the Hamiltonian. 
However the Hamiltonian approach necessitates a splitting of spacetime at least to some 
extent. One constructs a 3 -|- 1 foliation of spacetime by selecting a time function t such 
that the hypersurfaces, Sj, of constant t are space-like Cauchy surfaces. The standard 
Hamiltonian formulation for general relativity is the ADM representation (see. e.g., |1I|J23[ 



and PH] Ch 21), in which 4-covariant objects are decomposed into various 3-covariant parts. 



In particular, the spacetime metric, g^,y, is decomposed into the form, 

= -N'^dt^ + habidx'' + N''dt){dx'' + N^dt), (11) 

which depends on three spatially covariant parts: the lapse function N, the shift vector N°- 
and the spatial metric hab, induced on S^. The associated Hamiltonian density is obtained 
from 
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(12) 



where hat = dthat '■= £thab- Although this approach has led to much insight, it has the 
drawback that the resultant Hamiltonian formulation is only manifestly 3-dimensionally 
covariant. 

We prefer to use a more "covariant" approach to the Hamiltonian density. To this end 
it is convenient to use the so-called 'Palatini' method of treating the metric and connection 
independently. (Our approach is in many ways similar to that of Kijowski, see, e.g. 
Let us note some general features. First of all, since we are concerned with localization, 
we shall want to find the Hamiltonian which can evolve a finite spatial region. To achieve 
our 'covariant' formulation, we represent the time evolution direction as a covariant 4- vector 
field iV^. Quite generally the (4-covariant) Hamiltonian — which is essentially the Noether 
generator of translations (i.e.. Lie derivatives) along iV^ — is given by the spatial integral 
over (the finite spatial hypersurface) of a 4-covariant Hamiltonian density. In order 
to generate the Lie derivatives, the Hamiltonian density is necessarily linear in the time 
displacement vector field A^^ and its derivatives. Consequently it can be expanded into the 
form 

W{N) = N^'W, + d,[B^'''{N)l (13) 

where, it turns out that (at least for our representation) B^^{N) = — i3'^''(A^). 

On the other hand, beginning from the Lagrangian density, we can apply our Noether 
type argument to a translation along A^^, (see [0). Formally we then arrive at a conserved 
quantity (essentially the canonical energy-momentum density discussed earlier) which is 
actually just this same Hamiltonian density. From this analysis we learn a couple of impor- 
tant things. First, we find that ("on shell") the Hamiltonian density is necessarily conserved: 
dfjT-i^{N) = 0. At this point we want to call attention to the fact that the possibility of ad- 
justing the canonical energy-momentum density (^) by adding an automatically divergence 
free part exactly corresponds to adjusting B^'' in (|13]). Second, we find that H^^, is linear in 
the field equations and thus vanishes "on shell" . From this latter fact we conclude that the 
numerical value of the Hamiltonian is completely determined by the dB term, which, when 
integrated over a spatial region S^, via the divergence theorem, gives rise to an integral of 
B over the 2-dimensional spatial boundary dJ^t- The value of the Hamiltonian for a finite 
region is thus determined by the Hamiltonian boundary term and hence is quasilocal: it is 
associated with the closed spatial 2-surface boundary of the region. 

Now we turn to specific details for GR. Because we work with first order Lagrangians, we 
can easily obtain the Hamiltonian by merely rearranging the Lagrangian into the field theory 
analogue of L = p^q'" — H; essentially from (P we simply get C = P^dt4>^ — {—P%da(t)^ + A). 
We first consider the Hilbert Lagrangian which is given by the spatial integral of the Hilbert 
Lagrangian density (^). The spatial integrand can be expanded in the form 

CHN^^dT.^ = {r%.7r^^<5^:!; - H^(A^)} rfS^, (14) 

where dE^ := ^e^a/s-rdx^dx'^ dx'^ and our definition of r"/?^, which simply reduces to 
N^d^V^p^ = dtT^p^ in adapted coordinates, is given in general in appendix ^ along with 
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other details regarding our choice of representation. From after a straightforward cal- 
culation, without discarding any total derivative term, we obtained the explicit expressions 

= -Ir'-p^xtt^'C: - T'^puDySi^l (15) 

}3^h{N) = N^Tt^''T''f,^6Z- (16) 

Note that, as predicted, the (spatial hj^er-) surface density part Ti^i. vanishes 'on shell', 
since it is linear in the non-metricity, DaQ^v, and, with vanishing non-metricity, the cur- 
vature contractions reduce to the Einstein tensor, G^jy. Hence, as expected, the value of 
the Hamiltonian comes only from the boundary term which gives the quasilocal energy- 
momentum. 

The Hamiltonian symplectic structure can be found, as in the Lagrangian case, by varying 
the Hamiltonian (regarding it as a function of T"' and its conjugate momentum n"'^^ 5!:!^l^ = 
S^Ti^'^ — S^^Ti'^^^). This variation fits the general pattern 

Sn^{N)= (field equation terms) +dX^''{N). (17) 

The field equation terms include a set of initial value constraints and dynamical equations 
p5|,p]j23[| which may be used to calculate the evolution of the gravitational fields. Here our 
focus, however, is on the variational boundary term C'^'^ = —C^^ which reflects the symplectic 
structure — and the implicitly built in boundary conditions — of the physical system with 
respect to the particular Hamiltonian density under consideration. The variation of the 
spatial hypersurface part, N'^Hy in addition to the field equation terms, gives rise to the 
total divergence 



Combining this with the variation of the boundary term (piq), we find that, for the present 
(Hilbert) case the Hamiltonian variational boundary term, C, takes the explicit form 

C],'{N) = -2t:^^ 5V\p 5;[-iV''l, (19) 

showing that the 'control variable' is — similar to the Lagrangian case — (certain projected 
components of) the connection coefficients. 



Expanding out the Hamiltonian boundary expression (|lB), and using the Dg = field 



equation to express the connection coefficients in terms of the metric gives 

B^j,'{N) ^ K-'V^gN-^gP^'^T^^, ^ (2«:)-iv^iVVV^(5/3^?.7 - d^d^)- (20) 
This is in fact the superpotential which gives rise to the M0ller pseudotensor From 



this calculation we have acquired two insights: first, the M0ller pseudotensor is essentially 
a quasilocal object and, second, it really gives the energy — the value of the Hamiltonian 
— for the particular Hamiltonian which generates time displacements in the case that the 
connection is fixed on the boundary. 

However, there are some shortcomings in this Hamiltonian (aside from the obvious fact 
that the boundary term is not covariant, which we will remedy further below). First, al- 
though the boundary condition at infinity is simply F = 0, we must consider the rate of 
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approach. With the standard fall offs, in particular ST — > 0(l/r^), the boundary term in the 
variation of the Hilbert Hamiltonian will not automatically vanish asymptotically, indicating 
that the Hamiltonian is not different iable on the phase space of all asymptotically flat fields; 
hence one must actually give the explicit asymptotic functional form for F at each instant 
of time. Second, the actual value of the energy calculated from the boundary term for the 
Schwartzschild solution is not the expected value M but exactly half of that amount. 

This fact is closely connected with a well-known problematical feature of Komar's co- 
variant expression p3| . 



Bk'{N) := -V^D^^N^^ = -^{g^^l^d.N^^ + N''g^^-r^\,), (21) 

which is equivalent to the M0ller superpotential ( PP] ) when the components of iV^ are con- 
stant. Long ago it was noted that if the Komar expression is normalized to give the correct 
energy-momentum then it gives a value for the angular momentum which is twice that de- 
sired, conversely if it is normalized to give the correct angular momentum it then gives only 
half of desired amount for the energy. (The proper way to reconcile the M0ller-Komar su- 
perpotential with the desired energy-momentum and angular momentum results was found 



some time ago [^,^,0 . In fact the results of those works are forerunners of our preferred 
expressions discussed below.) Here we found that, from the standard normalization of the 
Hilbert Lagrangian, the associated boundary term in the Hamiltonian (obtained without 
discarding or modifying any boundary terms) naturally gives rise to the M0ller-Komar su- 
perpotential with the latter normalization. 



M0ller himself later noted that the Hilbert Lagrangian leads to his superpotential p? 



Also, it has long been known that the Hilbert Lagrangian leads via a Noether argument to 
the Komar superpotential (see, e.g., [Q). (The Komar potential has also been obtained in a 
Hamiltonian treatment via a Legendre transformation from the Einstein Hamiltonian [0.) 
Although the factor of 2 problem with Komar's expression has also long been known, yet it 
seems that only very recent works |2^JT^ have explicitly noted that the normalization arising 
directly from the Hilbert Lagrangian gives only half of the expected energy-momentum. 

There is a very simple cure to this problem of getting half of the desired value. Exploiting 
the freedom in the Hamiltonian that we noted above — the freedom to add a divergence 
free term to the canonical energy-momentum density without changing the fact that it is 
conserved — we can modify the Hamiltonian boundary term. Adjusting the Hamiltonian 
boundary term 'by hand' will not change the equations of motion but it will change the 
boundary conditions and the value of the quasilocal energy-momentum. Note that such an 
adjustment is essential if we wish to obtain a Hamiltonian which will be differentiable on 



the phase space of all asymptotically flat fields, as was nicely explained some time ago |^ 
in connection with the usual ADM formulation. (Indeed the usual approach is simply to 
discard all boundary terms on the way from the Lagrangian to the Hamiltonian and then to 
fix up the Hamiltonian boundary term in the end to produce the desired behavior.) Given 
this freedom, we could simply double the boundary term in the Hilbert Hamiltonian. This 
would certainly take care of the problem of getting only half the value for the energy, but 
the symplectic structure in the variation of the Hamiltonian would then be modified. The 
necessary boundary condition would then require the vanishing of 

{n^ {T''p,57rf'^5Z + 5T%,7r^'^5Z + 25r^^7r^'^5-|,-5^ ) } , (22) 
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which leads to a rather unattractive, comphcated boundary condition requiring the vanish- 
ing of a combination of 5tt and 5T and which, moreover, would not automatically vanish 
asymptotically — so the Hamiltonian would still not be differentiable on the space of all 
asymptotically flat fields. 

Let us now briefly consider the Einstein Lagrangian density (|^. We take the dynamical 
"coordinate" variable to be t^^"' . The associated Hamiltonian density can be found from 

CEN^'dT.^ = {vv^^r^.^^!; - n%{N)]dT.^, (23) 

(our general definition of tt is given in appendix 0; in adapted coordinates it simply reduces 
to dtir""^). The Hamiltonian density 1-Le still has the same ADM surface part Ti^p (but 
when varied it is now to be regarded as a function of tt""" and its conjugate momentum, 
^V'^^aS'^ + iP'^'^a^^ — r^^). However the Hamiltonian has now acquired a different boundary 
term given by 

B'-^iN) := N^T^^'^T^.dZ^r- (24) 

This Hamiltonian boundary term, which arose directly from the Einstein Lagrangian density 
without discarding or adjusting any exact differentials, is a familiar object. Using the metric 
compatible field equation to replace the connection by derivatives of the metric leads to a 
well-known form of the expression, 

B^E^N) = {2n)N^{-g)-^/^gxrd^{7T^^7i''^ - vr^^vr^^). (25) 



This is exactly the Freud superpotential |T9[ whose divergence gives rise to the Einstein 
pseudotensor. The spatial integral of the Einstein pseudotensor yields a value which is 
actually quasilocal, it is given by the integral of the Freud superpotential over the closed 
2-boundary of the spatial region. This is identically the same boundary integral and thus 
the same quasilocal value as is determined by the Einstein Hamiltonian via its boundary 
term. Extending the region to spatial infinity yields the total energy-momentum, now with 
the proper normalization. 

The boundary term in the variation of the Einstein Hamiltonian takes a form which 
differs from the Hilbert case: 



Cl%N) = 2^^^ 67c'^''6;];Np'^ . (26) 

From the symplectic structure of this Hamiltonian variation boundary term we learn that this 
Einstein choice corresponds to holding fixed the contravariant metric density. With the usual 
asymptotics this term will vanish at spatial (but not at future null) infinity, consequently the 
Einstein Hamiltonian is automatically differentiable on the phase space of all asymptotically 
flat fields (spatially, while at future null infinity one must specify the detailed functional 
asymptotic form of the metric to describe the radiation). The symplectic response, however, 
reveals a deficiency. Since it is some projected components of the connection, it is not really 
a covariant object. An improved result could be obtained from the Lagrangian density (^, 
but we have already seen the important ideas so, instead of elaborating that case, we will 
just go on to our final forms for the Hamiltonian boundary term in the next section. However 
before we do that let us make a few observations. 
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The role of the variational boundary term in the Lagrangian and Hamiltonian formu- 
lations are similar, in that in both instances one can adjust the boundary term to change 
the implicit boundary conditions. However in the Hamiltonian case there is an additional 
entirely independent and very strong motivation which draws our attention to the boundary 
term and moreover invites us to modify it. Because the Hamiltonian is conserved, its value 
has a physical significance not shared by the Lagrangian. This conservation property is 
preserved under modifications of the boundary term (equivalently, preserved under adding 
a divergence free term to the Hamiltonian density). In the case of the Hamiltonian for 
dynamic geometry, the entire value actually comes from the boundary term. 

One important way in which the boundary terms in the Lagrangian and Hamiltonian 
differ is that the former determines a boundary condition on the whole boundary of a space- 
time region whereas the latter determines only spatial boundary conditions. By adjusting 
both of them accordingly we can independently choose what is held fixed on the initial 
time hypersurface and at the spatial boundary (which in fact is convenient for the different 
types — Cauchy vs. Dirichlet/Neumann — of boundary conditions typically required on these 
surfaces) . This fact is related to another way in which the Hamiltonian boundary term issue 
differs from that of the Lagrangian. At the Hamiltonian level, we can make boundary terms 
which depend on the displacement vector field A^'^ in various ways. This allows for many 
more possibilities than those like (|^ and (|^) that are available at the purely Lagrangian 
level. Consequently there is a bigger need for a suitably restrictive criterion. 

The plain fact is that we can entirely ignore the boundary term which arises from the 
Lagrangian and simply change the Hamiltonian boundary term to anything we want. How- 
ever our choice is constrained if we wish to satisfy an important physical desiderata: namely 
to get the desired energy-momentum values for empty space, weak fields and at spatial and 
future null infinity. This requirement shows up only at the Hamiltonian level; it is easily 
dealt with at that level whereas in general it is not so readily satisfied by a judicious ad- 
justment of the boundary term back at the Lagrangian level. In fact this requirement forces 
us to adjust the Hamiltonian boundary term away from that naturally inherited from the 
Hilbert Lagrangian. Moreover, it actually fixes the form of the Hamiltonian boundary term 
— but only to linear order. Going beyond the linear order we can use our freedom to build 
in certain boundary conditions via the Hamiltonian variation symplectic structure. 

One consequence of this freedom is that, not only the superpotentials for the M0ller 
and Einstein pseudotensors, but in fact also the superpotentials associated with all of the 
other pseudotensors are likewise acceptable Hamiltonian boundary terms. Here we briefly 



outline the argument which we have presented in more detail elsewhere [T^. Consider 
the pseudotensor idea: a suitable superpotential H^"^ = H^^'^^^ is selected and used to 
split the Einstein tensor thereby defining the associated gravitational energy-momentum 
pseudotensor: 

1 



Ky/^mt,-" := -m^G," + -9A(Ar^if/"), (27) 

where we have inserted a vector field to make the calculation more nearly covariant. The 
usual formulation is recovered by taking the components of the vector field to be constant in 
the present reference frame. Einstein's equation, G^'^ = kT^'^, can now be rearranged into a 
form where the source is the total effective energy-momentum pseudotensor 
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(28) 



An immediate consequence of the antisymmetry of the superpotential is that T^'^ is a con- 
served current: d^[{—gy^'^Tf^'^] = 0, which integrates to give a conserved energy-momentum. 
The energy-momentum within a finite region 



-P(iV) := - / 



[iV^v^(-G/ - T/) 



1 



2k 



I B{N) ^ H{N), 

Js=dT, 



(29) 



is seen to be just the value of the Hamihonian. Note that Ti^fj, is the covariant form of the 
ADM Hamiltonian density, which has a vanishing numerical value, so that the value of the 
Hamiltonian is determined purely by the boundary term B{N) = —N^^{l/2K)H^'^^{l/2)dSy\. 
Thus for any pseudotensor the associated superpotential is naturally a Hamiltonian boundary 
term. Moreover the energy-momentum defined by such a pseudotensor does not really 
depend on the local value of the reference frame, it is actually quasilocal — it depends (through 
the superpotential) on the values of the reference frame (and the fields) only on the boundary 
of a region. 

The Hamiltonian approach endows these quasilocal values with a physical significance. 
To understand the physical meaning of the quasilocalization, calculate the boundary term 
in the Hamiltonian variation: 



TP) 



dS. 



TP- 



(30) 



(This result differs slightly from (|I^) because the ADM form of the Hamiltonian used here 
does not contain a term proportional to Dn.) For example for the Einstein pseudotensor, 
use the Freud superpotential (|20|) as the Hamiltonian boundary term in (p9l). Then the 
boundary term in the Hamiltonian variation has the integrand SItt^'^ N^^)r°' i^xS'^'^^, which 
shows not only that vr^"" is to be held fixed on the boundary, but also that the appropriate 
displacement vector field is N'^ = constant, and the reference configuration here is Minkowski 
space with a Cartesian reference frame. 

A minor variation on the preceding analysis results from choosing a superpotential with 
a contravariant index: H^'^'^ = }Jp-WM_ A further variation: H^'^" := df^H'^"'^/^, along 
with the symmetries H^'"'''^ = H''^^''^ = ifMI'^/^] and i^^'I"'^/^] = 0, leads to a symmetric 
pseudotensor — which then allows for a simple definition of angular momentum, see 
§20.2. We can cover these options simply by using the displacement vector field to make 
modifications like N^'H^''^ — > Nf.H^"^. 

In this way we see that each of the pseudotensors actually gives the value of the quasilocal 
energy-momentum for an acceptable Hamiltonian. In each case, via the Hamiltonian bound- 
ary variation symplectic structure, this quasilocal energy-momentum is associated with some 
definite physical boundary conditions |T0|. Note that this same type of argument extends to 
superpotentials (i.e., Hamiltonian boundary terms) that are more general than the classic 
linear-in-displacement form associated with the traditional pseudotensors. In particular one 
can include first (and even higher) derivatives of the displacement, as occurs in the Komar 
expression 
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In summary, similar to the Lagrangian analysis which we have discussed, the boundary 
term in the Hamiltonian variation, C, generally does not vanish, so the Hamiltonian is not 
differentiable for general field values. A modification, achieved by adding a total derivative 
term to the Hamiltonian, adjusts B without changing the field equations and can compen- 
sate, making C vanish for suitable preselected boundary values. The exact form of such 
an adjustment still has infinite possibilities, this allows for an infinite number of different 
gravitational energy definitions However, each of them has its own unique expression 
for the Hamiltonian boundary variation C. The symplectic structure of this term reveals the 
implicit boundary conditions and thereby gives a physical interpretation for each quasilocal 
energy-momentum expression. Thus, for each well-defined Hamiltonian boundary expres- 
sion, one can, via the Hamiltonian analysis, find its associated symplectic structure which 
shows the built in control mode, or equivalently, the implicit boundary conditions. 

III. QUASILOCAL ENERGY-MOMENTUM 

Here we describe our Hamiltonian boundary term expressions for quasilocal energy- 
momentum. Our major tool is the symplectic analysis of the Hamiltonian boundary vari- 
ational principle. We associate each possible Hamiltonian boundary term expression with 
the boundary conditions identified via the symplectic structure of the boundary term in 
the variation of the Hamiltonian. There are an infinite number of possible Hamiltonian 
boundary terms, and correspondingly an infinite number of possible boundary conditions. 
We greatly reduce this infinity by applying a covariance criteria. 

In the previous section we saw that the Hilbert Hamiltonian had problems asymptoti- 
cally while the Einstein Hamiltonian gave good asymptotic values but had a non-covariant 
response, being linear in the connection. These shortcomings necessitate, as we saw at the 
Lagrangian level, the introduction of a reference geometry. Hence for regulating the varia- 
tional boundary term, a background manifold with a suitable geometry, {M,g^i,,T°'^i,), is 
introduced as a reference configuration. The gravitational energy-momentum is understood 
to be measured with respect to this selected background. Any modification of the Lagrangian 
or Hamiltonian boundary term changes the symplectic structure and the boundary condi- 
tions. Here, from the two examples we considered, we obtain modified versions of Bh and 
Be which have the same the control modes, respectively T^ap or tt^'^ , but their responses 
are given an improved "covariant" form. 

For the metric density (in deference to the traditional choice of variables we refer to it as 
the "Dirichlet") control mode, the background is just what we need to make the responses 
become tensorial objects without changing the control variables. Its symplectic structure in 
the variational boundary term is required to have the form 

C:'{N) = 2Ar%^ ^vr^^^^^t^AT"], (31) 

where the A means the difference of variables between physical and background configura- 
tions (i.e., AT^ajS '■= — and Avr^*^ := ix^''^ — ti^'^). Now the response is a combination 
of AF which is a tensor. Moreover the whole Hamiltonian boundary variation term is now 
the projection along the displacement vector field of a four dimensionally covariant object, 
a vector density which vanishes asymptotically (spatially) with standard fall offs — showing 
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that the Hamiltonian is differentiable on the space of all asymptotically flat fields. In order 
to obtain this desired C^, we must modify Be- The modified quasilocal energy-momentum 
boundary term, S^, was found in [O-M to be 



B^%N) = N^nf^'AT^^SZl + N^T'^ ^rAn^^SZ- (32) 

Similarly, for our connection ("Neumann") control mode, the boundary term in the 
Hamiltonian variation can also be improved by incorporating reference quantities in the 
form 

Cp^ = -2A7f^'' Sr^'^f^S^j^NP^ . (33) 

This symplectic expression is again the projection along the displacement vector field of a 
four-covariant vector density which automatically vanishes asymptotically (spatially, with 
standard fall offs) indicating that the Hamiltonian is differentiable on the space of all asymp- 
totically fiat fields. This version follows from the adjusted Hamiltonian boundary term 

= N^nf^^AVp^d^ll + N''T''f3^A7r'^^5^l (34) 

Note that the two modes are complimentary: the Hamiltonian boundary variation symplectic 
relation for one can be obtained from the other just by interchanging the control- response 
roles. 

From the variables at hand there are two other Hamiltonian boundary term expressions 
which can be constructed: 

= N^TT^^AT'^f.^SZl + N^T'^f.^Air^^SZ, (35) 
= N^7rf^^AT''p^5^ll + N^T" ^^Anf^^S^l (36) 

Bq has the interesting property of being linear in the dynamic variables vr, F while Bi is 
linear in vf, F. The two associated Hamiltonian variation boundary terms (both of which 
automatically vanish asymptotically with standard spatial fall offs) have a remarkable A ^ 6 
symmetry: 

C^o' = -ST"p,A7i''''m6l^^ - AT'^p,N^67i''''6Z, (37) 
CI" = AF"^^57r'5'^iV'^5^e]; + 6r%^N^A7r^^6Z. (38) 

However, neither has the J^'^N^'^ form of a projection along A^^ of a 4-dimensionally covariant 
vector density, only our two expressions (^,^) (or constant linear combinations thereof) 



leading to (31,33) have this desirable 'covariant' property. 

Returning to our two 'covariant' expressions, there is a technical hitch here that needs 
discussion. Although the Hamiltonian variation control-response symplectic structure has 
a nice covariant form, the Hamiltonian boundary terms themselves ( p2| , p4D are not fully 
covariant. This is an inevitable consequence of our particular style of first order 'independent 
metric, frame and connection' formulation, as we briefly explain here (the main technical 
point is that we actually treat the connection as a one form; for further remarks see appendix 
The connection is not a covariant object. The Hamiltonian must generate the evolution 
of the connection coefficients including the reference frame gauge dependent part (which 
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depends on the displacement vector field differentially). This latter task is the duty of 
the NTDir term in the Hamiltonian. The Hamiltonian boundary term then includes an 
associated piece with the form A^FAtt. The contribution of this piece to the value of the 
energy-momentum is a mixture of a covariant physical contribution along with an energy- 
momentum associated with the particular reference frame. Fortunately these contributions 
can easily be separated by using the identity 

N^T^p^ = Df,N" - dpN^, (39) 

to replace the A^F terms. The ON terms produce a noncovariant (reference frame dependent) 
unphysical contribution (which can usually be made to vanish in a specially selected frame) 
and should be dropped (for the purposes of calculating physical energy-momentum but not 
for calculating the evolution equations). This leads to the final fully covariant form of our 
Hamiltonian boundary quasilocal energy-momentum expressions: 

(iV) = iV^vr^^AF^^^CAr + ^/siV-Avr'^^CA, (40) 
Bf^'iN) = N^Tr^^AT'^p^d^l'l + Df^N'^Ajr^^SZ- (41) 

We wish to emphasize that an alternate, fully covariant, c?zrect derivation of these expressions 
can be obtained from a different representation as indicated in Appendix 0. 

After a bi- metric manipulation (see Appendix^), the above expressions can be rewritten 
in the following compact and remarkably similar forms: 

i3f = 2A{tt^^''DxN^'^) + N^k^'in) - A^^r(7r), (42) 
= 2A{7T^^''DxN^'^) + N^k^'in) - A^^r(7r), (43) 



where 



F(7r) := TT^'^AF^, - tt^^AF^/^, (44) 



and ^'^(vr) has the same form with vr replacing tt. 

At first we were surprised to learn that our Dirichlet expression (HO) is exactly identical 



with an expression obtained by Katz et al. p^,p8|,p9| which was derived in a completely 



different way, namely by applying the Noether conservation theorem to the Lagrangian 
density (compare with (^): 

= -7r^'^(AF\,AF%, - AF\,AF^,) + Att'^'^R^,, (45) 

which includes background terms in addition to terms quadratic in the first derivatives of (^^j,. 
In retrospect we realize that our having found an identical energy-momentum expression is 
not so surprising after all. The Hamiltonian approach, as we discussed, is closely connected 
with the Noether approach. Moreover our covariance requirement leaves little room in the 
Hamiltonian boundary term for anything else except expressions that can be inherited from 
a suitable four dimensionally covariant Lagrangian. 

Comparing the remarkable similarity in the form of our alternate Neumann expression 
(PI) with that of (0), invites us to consider also obtaining it from a Lagrangian density. 
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The desired result is obtained simply by interchanging the roles of g and g (consequently 
A — i> —A) followed by an overall sign change. Thus, we found that (^) can be derived 
by the same Noether argument used by Katz and coworkers from the following Lagrangian, 
which is quadratic in the first derivatives of g^i,: 

= -TT^^AT^xAT^^, - AT\pAT^x) + An^'^R^,. (46) 

Via the Hamiltonian boundary term symplectic structure we identified the two quasilocal 
energy expressions (5^,il) as corresponding to Dirichlet or Neumann type boundary con- 



ditions respectively. Comparing their respective Lagrangians we see the relation between 
these two expressions from a new point of view. There is an amazing symmetry relating the 
two expressions: one passes into the other simply by interchanging the role of the dynamic 
physical and background variables. Hence we can regard our "Neumann" expression ( ^Tj ) 
as giving the energy-momentum of the "reference" geometry measured with respect to the 
"dynamic" geometry using "Dirichlet" boundary conditions. And, likewise, our "Dirichlet" 
expression gives the "Neumann" energy-momentum for the "reference" geometry compared 
to the "dynamic" geometry. This symmetry, however, has an intriguing asymmetry: one 
might have expected the energy-momentum of the dynamic space referenced to the back- 
ground to have the same magnitude when the roles are reversed without any reversal in 
the type of boundary condition (on the other hand, one could argue that asymmetries are 
common when the reference point is changed, e.g., going from 4 to 5 and back to 4 can be 
described as a 25% increase followed by a 20% decrease). We suspect that there is some, as 
yet unidentified, underlying principle which could have been used to anticipate this curious 
symmetry and asymmetry. 

Our two boundary expressions are not the only ones for gravitational energy-momentum. 
They are simply the only ones which satisfy our covariant Hamiltonian symplectic structure 
criterion. Covariance is a very important property; we believe that a covariant theory should 
have covariant quasilocal energy-momentum. Nevertheless it is well to keep in mind that 
some other property may be regarded as even more desirable (also, perhaps our particular 
implementation of the covariance requirement could be generalized). Then one could ex- 
ploit the freedom in selecting the Hamiltonian boundary term to achieve a different goal. For 
example Kijowski and Jezierski [^,0 have used the constraints and the actual boundary 
conditions required by the field equations to identify and control certain variables repre- 
senting the true physical degrees of freedom. This necessitates decomposing the dynamic 
fields into various space, time and boundary components. Consequently their expression 
for quasilocal energy-momentum is not covariant (although it would be interesting to try to 
recast it into that form). More recently, following along the lines of Rosenfeld and Belin- 
fante, Petrov and Katz have used what amounts to the same Hamiltonian boundary term 
freedom that we have exploited to achieve a "symmetric" energy-momentum expression |^ 



To any such alternate expressions one can apply, just as we did for the pseudotensors, our 
Hamiltonian boundary variation symplectic analysis to reveal the implicit spatial boundary 
conditions. 



16 



IV. CONCLUSIONS 



In summary, variational principles can yield not only field equations but also boundary 
conditions; the latter can be modified by adding a total derivative, which is equivalent to a 
boundary term. The Hamiltonian boundary term for dynamic spacetime governs the value 
of the Hamiltonian. For each finite region its value yields a quasilocal energy-momentum. 
The boundary term in the variation of the Hamiltonian has a symplectic structure, which is 
uniquely determined by the choice of quasilocal expression. Requiring it to vanish associates 
to each different quasilocal expression distinct boundary conditions. This approach provides 
a physical interpretation for many of the well-known gravitational energy- momentum expres- 
sions including all of the pseudotensors, associating each with unique boundary conditions. 
Among the infinite possibilities, we found only two Hamiltonian-boundary-tcrm quasilocal 
expressions which correspond to covariant boundary conditions; they are respectively of 
the Dirichlet or Neumann type. Our Dirichlet expression coincides with the expression re- 
cently obtained by Katz and coworkers using Noether arguments and a fixed background. 
A modification of their argument yields our Neumann expression. 

Some key points we have noted in our analysis are: 

• The boundary-variation-symplectic-structure principle connects the choice of bound- 
ary term with boundary conditions. The Lagrangian boundary term can be adjusted 
to affect a canonical transformation. It governs the boundary conditions on the 3- 
dimensional boundary of a spacetime region, including the initial time spacelike hy- 
persurface. 

• The Hamiltonian boundary term governs the boundary conditions on the 2-dimensional 
boundary of the spatial region at each instant of time. The value of the Hamiltonian 
for dynamic geometry theories including general relativity is determined entirely by the 
Hamiltonian boundary term. It gives the quasilocal energy-momentum. Our freedom 
to adjust the Hamiltonian boundary term is justified by the conservation law. The 
Hamiltonian boundary term depends on the displacement vector, hence it has (in 
principle) more freedom than is available at the Lagrangian level. However the value 
of the Hamiltonian, the energy-momentum, also has physical 'correspondence limit' 
constraints which have no analog for the Lagrangian. The boundary term freedom 
we exploit here is essentially the same freedom used in constructing 'new improved 
symmetric energy- momentum tensors'. 

• The Einstein and M0ller (Komar) pseudotensors arise quite naturally, but the latter 
has several more shortcomings. All of the pseudotensor superpotentials are possible 
Hamiltonian boundary terms. Consequently all pseudotensors have quasilocal energy- 
momentum which is identical to the value of the Hamiltonian for an acceptable choice 
of boundary term, which, in turn, corresponds to some definite boundary conditions. 

• Our 'covariance' criterion removes most of the freedom (leaving only two choices). In 
hindsight we sec that it essentially restricts us to expressions which could be obtained 
(without any adjustments by hand) by projecting a judicious choice of Lagrangian 
boundary term. Our Dirichlet mode unexpectedly coincides with that of Katz et al. 
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In retrospect that is no surprise. Our Neumann mode can be interpreted as the Katz 
et al. energy-momentum of the reference geometry referred to the dynamic geometry. 



We note some features of our quasilocal energy-momentum expressions: 

• We found only two "covariant" Hamiltonian boundary expressions. They each give 
rise to a boundary term in the variation of the Hamiltonian which has the form of 
a projection of a covariant vector density along the displacement vector field. The 
form of this variation boundary term shows that the respective Hamiltonians evolve 
field values with Dirichlet or Neumann type boundary conditions. With standard fall 
offs, the two Hamiltonians have well defined variational derivatives on the space of 
asymptotically flat fields at spatial infinity. 

• Our expressions depend on a reference configuration, which is required only on the 
boundary. The reference configuration determines the zero point for all of the quasilo- 
cal quantities. The obvious choice is Minkowski space; alternatives which may be 
more appropriate for certain applications include (anti-)de Sitter space, a Friedmann- 
Robertson- Walker cosmology, and Schwartzschild geometry. Some options for attach- 
ing an appropriate reference configuration to a dynamic boundary were discussed in 
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Our expressions also depend on a displacement vector field which selects the associated 
component of the quasilocal energy-momentum. How to choose the exact form of this 



vector field was discussed in the recommended choice is a Killing vector of the 
reference geometry. In addition to energy-momentum (obtained from a spacetime 
translation), for a suitable choice of rotational displacement, the expressions also give 
angular momentum. 

Our expressions reduce to expressions proposed by others in the appropriate limits, in 
particular to the well known quasilocal expressions of Brown & York |P and asymp- 
totically to that of Beig & 6 Murchadha 0. Asymptotically they are equivalent to 
an expression which gives the expected values at spatial infinity (for asymptotically 
fiat and ajoti-de Sitter solutions) [^]. Moreover, asymptotically, at future null infinity, 
our Dirichlet expression yields the expected Bondi values p2|. Quasilocally, we have 
evaluated them for spherically symmetric spacetimes [|1^,0. 

Katz and coworkers have applied their expression (equivalent to our Dirichlet expres- 
sion) at future null infinity to cosmology [^] and Mach's principle We have 



applied our formulation to black hole thermodynamics ||T3|JT^ to obtain the first law 
and an expression for the entropy. 



More generally our work reveals some of the merits of the symplectic Hamiltonian bound- 
ary variational principle. In particular it allows us to supplement the usual (correspondence 
limit to weak field and asymptotic forms) constraints on quasilocal energy-momentum ex- 
pressions with a principle which connects each quasilocal expression with a distinct boundary 
condition. Coupled with the covariance criteria the form of the quasilocal energy-momentum 
expression is then strongly restricted. 
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APPENDIX A: DYNAMICAL DETAILS 

Our Hamiltonian formalism is adapted to evolving the components of objects including 
the connection coefficients, hence it includes some (unphysical) dynamic reference frame 



gauge generation features. There are alternate representations (e.g., ||3T|) in which these 
terms do not show up. 

Our general approach is to work with a dynamical Lagrangian and Hamiltonian formu- 
lation which gives independent equations for evolving the frame, metric and connection, 
and which handles a wide range of theories, including geometric gravity theories and gauge 



theories in a uniform way [|14[. Note that in such a formalism the Hamiltonian must include 
the ability to generate a general time dependent (purely gauge) evolution of the frame and 
the associated induced effects on the components of geometric objects. In general we found 
that there are certain technical advantages in using a differential form representation. How- 
ever in the present work we wanted to make our results for the specific case of Einstein's 
GR more accessible to others, so we transcribed it into the ordinary "holonomic frame" 
representation. To achieve this some choices must be made: in particular, how to deal with 
the metric, frame and connection variables and how to impose the vanishing torsion and 
metric compatible constraints. We want to keep our first order form, so we certainly need 
the connection and metric to be independent at least to some extent. We elected to impose 
vanishing torsion 'a priori' and thus to use a symmetric connection and a variational princi- 
ple which would give the metric compatible condition as a (vacuum) field equation. Because 
we are using a holonomic frame the evolution of the frame is rather trivial, so we dropped 
it and its conjugate momentum from our list of dynamic variables. Nevertheless we did not 
want to depart far from the form of our earlier more general work. Thus the expressions 
given here are actually obtained by specializing our earlier work. In particular our time 
derivative is specified by projecting the Lie derivative £n '■= diN + iNd of components of the 
connection one-form and its conjugate momentum 2- form: 

t'^fs^dx^ := £Nir''pxdx^) = (Nf'd^n^y + np^d^N^')dx\ (Al) 
TT^'^e^, := i?^(7r^'^e,.) = (iY^^Avr^'^^^ + ^t:^^ d.N'^dZDep^, (A2) 

where e^y := {l/2)t^yai3dx'^ A dx^ . These "definitions" differ from the usual holonomic 
expression of the components of the Lie derivative for the contravariant metric density: 

^jvvr'^^ := dxiN^T^^""") - Tr^^^S^iV^ - vr^^S^iV^ (A3) 

and the connection coefficients: 

£^rV := + D^DpN'^. (A4) 
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The difference, however, shows up only in terms proportional to the derivative of N. (All 
such terms vanish if the coordinates are adapted so that £iq reduces to dt). A feature of our 
approach is a frame-gauge generating term in the Hamiltonian density and an associated 
term in the Hamiltonian boundary quasilocal expression. We like this set up for it is just 
the way things come out for the vector potential (i.e. connection one-form) in gauge theories 
like electromagnetism and Yang-Mills. 

It is certainly possible to use the usual Lie derivative in a Hamiltonian formulation 
(see, e.g., ||T6| , |32|j3T| ) . We choose to avoid it also because it includes an inconvenient for 
us second derivative of (which necessitates adjustments in our argument regarding the 
form of and the vanishing of the Hamiltonian density). A price we pay is then an awk- 
ward term like N'^T°' pr^'^^'^S'^x in each of our quasilocal energy-momentum Hamiltonian 
boundary expressions. We have argued that such terms are necessary to give us the Hamil- 
tonian evolution and boundary variation symplectic structure in our representation. Un- 
fortunately the quasilocal energy-momentum, defined as the value of our Hamiltonian, con- 
sequently includes both a physical and an unphysical, reference frame gauge dependent, 
contribution. To separate these effects we rearrange the symmetric connection identity 
{£NepY = [N, epf = (Vjve/3 - V/^iV)" to give 

Ar'^rV = D^Ar" + (i:jve/3)", (A5) 

which can be used to replace the A^F factors. The £NGf3 term is a non-covariant, dynamic 
reference frame piece. Its contribution to the energy-momentum can be thought of as an 
energy associated with the observer. In fact, for any given displacement vector field A", we 
can choose the reference frame e^j so that it vanishes. 

Having introduced this identity, an alternative approach is available. We could treat this 
term in the same way as its analogue is dealt with in other representations, in particular 
Kijowski's [pri. Note that, since it includes a time derivative, it really has no place in a 



Hamiltonian. Rather it should be treated as term belonging to the p^q'' part of the action, 
a term that shows up in a 2-dimensional integral over the boundary of the spacelike hyper- 
surface rather than in the 3-dimensional hypersurface integral. An easy way to establish 
the association between this part of our representation and Kijowski's is to consider the 
frame to be orthonormal. Then its time evolution is just an instantaneous Lorentz boost (in 
the spacetime 2-plane orthogonal to the spatial boundary) by a hyperbolic angle a6t. The 
associated 'conjugate momentum' is the area of the 2-surface. Hayward |^ gives another 



route to time derivative terms on the spatial 2-boundary. He uses the fact that the total 
boundary term in the Einstein action (|^) can be expressed as the extrinsic curvature of the 
boundary. The standard definition of the extrinsic curvature involves the normal to the 
boundary surface. But converting the total derivative form to a surface integral is then a 
delicate task, as the normal is discontinuous on the corners of the usual 3-boundary, which 
consists of an initial and final constant time spacelike hypersurface connected by a topolog- 
ically S*^ x [ti,tf] type 3-manifold. This leads to contributions in the action given by the 
difference between an integral over the final and initial 2-boundary. Contributions which 
can, in turn, be written as the integral over time of a total time derivative of a 2-boundary 
term. 

Actually it is not difficult to obtain a fully covariant Hamiltonian density with our 
fully covariant quasilocal boundary terms. Beginning from the Hilbert Lagrangian (13), the 
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Hamiltonian can be derived by using the usual Lie derivative of a connection 



W{N) ■= £Nr%y''5Z - N^'Ch (A6) 
= -liV-i?°^^,7r'^^5^^^J - DpN'^dZDuT^^'' + (DpN'^Ti^'^dZ) • (A7) 

The boundary term here is just the Komar superpotential (with the normahzation that 
gives half of the desired energy- momentum) . The boundary term in the variation of this 
Hamiltonian will not automatically vanish asymptotically; hence this Hamiltonian requires 
the explicit functional form of the connection to be fixed on the boundary even asymptoti- 
cally. Consequently this Hamiltonian should be adjusted. Replacing the boundary term by 
one of the improved boundary terms (^OD or (|4l|) gives a fully 4-covariant Hamiltonian for 
general relativity. Explicitly calculating the resultant boundary term in the variation of the 
Hamiltonian then leads to the desirable asymptotically well behaved covariant symplectic 
structures (|31|), (p^). For constant components A^'^ these fully covariant Hamiltonian density 
plus boundary term expressions reduce to (OJl^JS^JSl) . 



APPENDIX B: GEOMETRY OF BI-METRIC SPACETIME 

A background is needed to determine well-defined conserved quantities in GR. For the 
special choice of mapping and coordinates such that a point P of the physical configuration 
is mapped into a point P of the background and both are given the same coordinates x^, 
the whole system can be looked at as a spacetime M possessing two metrics g^j_y and ^^i,. 
Geometric quantities with respect to each metric can then be reformulated in terms of the 
difference between them. In particular each metric determines its own associated Levi-Civita 
connection and Riemannian geometry. 

The simplest case for the connection, from which all others can be derived is 

(d^ - D,) = Ar\,N\ (Bl) 

where the variables and operators are denoted with or without a bar consistently with the 
notation for the metric, and the symbol A means the difference of operands between two 
metrics such as AF = F — F. This identity shows that AF, being the difference between two 
connections is a covariant tensorial object. 

The Ricci tensor R^^i, {RpLv) with respect to F^^ (f can be rewritten, respectively as 

R^, = 5aAF\, - D,.AF\a + AF''^,AF\a - AF^^aAF^p + (B2) 
R^, = -Z^aAF\, + D^AF\a + AF^^,AF\a - AF^^a AF^^ + R^,. (B3) 

Two other useful identities concern the total derivative terms, which are added to the 
Hilbert Lagrangian density in order to make the Lagrangian density quadratic in the first 
order derivatives of the metric: 

a^F(7r) = -7r'^'^{(DAAF\, - D^AF\a) + 2(AFVAFV - AF%aAF\,)}, (B4) 
d,k^{7,) = -7r^''(DAAFV - ^M^r^A), (B5) 



where the k^^iji) and k^(TT) were defined in connection with ( ^41) 
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